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ABSTRACT 


Measurements  of  atmospheric ' transmission  during  periods  of  snow 
fall  were  made  in  the  visible  and/4nf raredSduring  March  and  April  1985 
The  measurements  are  compared  to  transmission  predicted  by  the  XSCALE 
module  of  the  Electro-Optical  Systems  Atmospheric  Effects  Library 
(EOSAEL).  'The  XSCALE  predictions  correspond  well  with  measured  trans- 
mi  ttances  in  the  visible  but  predictions  in  the  IR  are  considerably 
less  than  the  measured  values*  This  suggests  that  the  snow  transmit¬ 
tance  algorithm  used  in  XSCALE  is  more  accurate  for  transmissometers 

whose  receiver  diameters  are  much  larger  than  the  beam  diameter. 

_  .  ...  f  - 


RESUME 


^  r  - 


La  mesure  de  la  transmission  atmospherique  dans  le  visible  et 
1' infrarouge  a  ete  effectuee  au  cours  de  chutes  de  neige  pendant  la 
periode  mars-avril  1985.  Les  mesures  sont  comparSes  aux  transmissions 
predites  par  le  module  XSCALE  du  "Electro-Optical  Systems  Atmospheric 
Effects  Library"  (EOSAEL).  Les  predictions  de  XSCALE  correspondent 
bien  aux  transmittances  mesurees  dans  le  visible  mais  les  predictions 
dans  l'IR  sont  considerablement  moindres  que  les  valeurs  mesurees. 

Cela  suggere  que  l'algorithme  utilise  pour  calculer  la  transmittance 
dans  la  neige  est  plus  precis  pour  les  transmissom^tres  dont  le  recep- 
teur  est  beaucoup  plus  large  que  le  diametre  du  faisceau. 
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LIST  OF  SYMBOLS 


DESCRIPTION 


Fraction  of  total  extinction  due  to  snow 

Length  of  transmission  range 

Number  density  of  scattering  particles 

Scattering  efficiency  of  a  particle 

Perceived  reduction  in  Q  due  to  forward  scattering 

Effective  extinction  efficiency 

Effective  extinction  efficiency  averaged  over 

transmission  range 

Relative  humidity 

Temperature 

Visibility  or  meteorological  range 

Distance  from  receiver  along  transmission  path 

Radius  of  scattering  particle 

Mean  snow  particle  radius 

Radius  of  transmissometer  receiver 

Volume  extinction  coefficient 

Effective  extinction  coefficient  along  transmission 
path 

Fog  extinction  coefficient 
Snow  extinction  coefficient 
Extinction  coefficient  In  the  visible 
Particle  size  parameter 
Angle  subtended  by  receiver 
Wavelength  of  light 

Phase  function  of  scattering  particles 
Total  atmospheric  transmittance 
Transmittance  due  to  atmospheric  aerosols 
Transmittance  in  visible  band 
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1.0  INTRODUCTION 


The  Electro-Optical  Systems  Atmospheric  Effects  Library  (EOSAEL) 
is  a  library  of  computer  codes  developed  at  the  US  Army  Atmospheric 
Sciences  Laboratory  (ASL)  which  predict  the  effect  of  natural  and  bat¬ 
tlefield  gases,  aerosols,  smokes  and  dust  on  electro-optical  system 
performance.  EOSAEL,  which  is  divided  into  20  subprograms  or  modules, 
calculates  extinction  coefficients,  atmospheric  effects  on  radiation 
propagation  and  contrast  reduction  for  wavelengths  in  the  visible, 
infrared  and  millimeter  wave  regions.  XSCALE  is  the  EOSAEL  module 
which  calculates  visible  and  IR  transmittance  under  low  visibility 
weather  conditions,  such  as  fog,  rain  or  snow. 

In  order  to  estimate  transmittance  through  falling  snow,  the 
XSCALE  algorithm  requires  as  input  the  temperature,  relative  humidity, 
receiver  radius  and  the  visibility.  An  estimate  of  the  mean  snow  par¬ 
ticle  radius  is  made  based  on  the  temperature,  and  the  visibility  is 
used  to  calculate  the  extinction  coefficient  along  the  transmission 
path.  Since  snow  particles  scatter  strongly  in  the  forward  direction, 
the  measured  transmittance  depends  not  only  on  the  extinction  coeffi¬ 
cient  but  also  on  the  size  of  the  transmissometer  receiver.  The  XSCALE 
algorithm  is  designed  to  take  the  effect  of  forward  scattering  into 
account  based  on  the  receiver  radius. 

In  this  report,  estimates  of  transmittance  in  the  visible  and  in 
the  3-5  um  and  8-12  ym  bands  are  made  using  XSCALE  and  compared  to 
actual  transmittance  measurements  made  at  DREV  under  snowfall 
conditions . 

This  work  was  performed  between  March  1985  and  January  1986 
under  the  terras  of  a  Memorandum  of  Understanding  (MOU)  between  ASL  and 
DREV.  The  work  was  performed  under  PCN  33A11,  Aerosol  Studies. 
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2.0  INSTRUMENTATION 

2 . 1  Transmission  Range  and  Locale 

The  transmission  range  is  located  at  Defence  Research  Establish¬ 
ment  Valcartler  (DREV),  situated  approximately  20  km  northwest  of 
Quebec  City.  The  local  environment  is  rural.  The  receiver  station  is 
located  on  the  top  floor  of  a  laboratory  building  with  the  receivers 
9.4  m  above  the  ground.  The  sources  are  sheltered  in  a  hut  at  ground 
level,  538.4  m  from  the  receivers  along  the  13°  NW  direction.  The 
transmissoraeter  path  is  slightly  slanted  at  an  angle  of  1.43°.  During 
the  measurement  period  (between  January  and  April  1985),  the  land  below 
the  transmission  path  was  completely  snow  covered.  The  meteorological 
and  aerosol  sensors  were  located  on  the  roof  of  the  receiver  station 
approximately  3  m  above  the  receivers. 

2 .2  Visible  Transmissometer 

The  visible  transmissometer  consists  of  a  10-mW  HeNe  laser 
source  and  a  Fresnel  lens/silicon  detector  receiver.  The  laser  beam  is 
chopped  at  a  frequency  of  480  Hz  and  collimated  to  a  5-cm-diaraeter 
beam.  The  demodulation  reference  signal  is  derived  from  the  chopper 
and  is  transmitted  to  the  receiver  station  by  a  small  FM  radio 
transmitter. 

The  receiver  consists  of  a  25.4-cm-diaraeter  f/0.8  Fresnel  lens, 
a  focal  plane  aperture  and  a  5  x  5  cm2  silicon  detector  positioned  1  cm 
behind  the  aperture  (see  Fig.  1).  The  focal  plane  aperture  is  used  to 
limit  the  field  of  view  of  the  receiver  to  40  mrad.  The  detector  out¬ 
put  is  connected  to  a  lew  input  impedance  current-to-voltage  converter 
which  provides  a  voltage  output  proportional  to  the  intensity  of  the 
incident  radiation.  Tests  have  shown  that  the  receiver  response  to  the 
modulated  laser  light  is  Independent  of  the  level  of  ambient  light. 
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2.3  Infrared  Transmissometer 


The  source  of  the  DREV  IR  transmissometer  is  a  Barnes  model 
14-708  IR  source  consisting  of  a  650°C  blackbody,  90-Hz  chopper  wheel 
and  12.2-cm-diameter  off-axis  parabolic  mirror  which  collimates  the  IR 
beam.  The  demodulation  reference  signal  is  derived  from  the  chopper 
wheel  and  is  transmitted  to  the  receiver  station  by  radio. 

The  receiver  was  developed  at  DREV  and  is  shown  in  Fig.  2.  The 
detector  is  an  Infrared  Industries  InSb/HgCdTe  sandwich  detector.  The 
radiation  is  collected  by  a  9.14-cm-diameter  off-axis  parabolic  mirror 
aid  focused  to  a  point  in  front  of  a  3-mm-diameter  germanium  field 
lens.  The  receiver  field  of  view  (FOV)  is  6.0  mrad.  The  lens  has  an 
antireflection  coating  for  the  8-14  pm  band  to  enhance  transmission  to 
the  less  sensitive  HgCdTe  detector.  The  IR  beam  is  intercepted  near 
the  focus  by  a  filter  wheel  containing  five  filters  of  bandwidths: 

3. 4-4.1  pm,  4. 5-5.0  pm,  3. 0-5.0  pm,  8-12  pm  and  10-12  pm.  For  this 
report,  the  measurements  made  with  the  relatively  wide  3. 0-5.0  pm  and 
8-12  pm  filters  will  be  used.  The  filters  are  advanced  sequentially 
under  the  control  of  a  microcomputer  which  also  controls  a  data  acqui¬ 
sition  system.  A  transmission  measurement  is  made  using  each  filter 
approximately  every  50  s.  The  calibration  of  both  transmissometers  is 
described  in  Ref.  1  and  the  absolute  accuracy  has  been  calculated  to  be 
3%  for  the  IR  transmissometer  and  1.5%  for  the  visible  transmissometer. 
The  main  specifications  for  the  transmissometers  are  summarized  in 
Table  I. 


2.4  Meteorological  Parameters 

Temperature,  dew  point,  atmospheric  pressure,  wind  speed  and 
wind  direction  were  measured  using  standard  meteorological  instruments. 
The  sensors  were  mounted  on  the  roof  of  the  receiver  station  and  the 
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TABLE  I 


Specifications  of  IR  and  visible  transmlssometers 


Transmlssometers 

Quantity 

Units 

Visible 

Infrared 

Collimator  objective  dia. 

cm 

5  ±  .1 

12.2  ±  . 

Receiver  objective  dia. 

era 

25.4  ±  .1 

9.14  ±  . 

Collimator  focal  length 

cm 

20  ±  .1 

65  ±  . 

Receiver  focal  length 

cm 

20  ±  .5 

63.4  ±  . 

Source  $  angle  FOV 

mrad 

.05  ±  .01 

5  ±  . 

Receiver  §  angle  FOV 

mrad 

20  ±  1 

3  ±  . 

Range 

m 

538  ±  1 

538  t  1 

Absolute  accuracy 

% 

1.5 

3 

Some  inconsistencies  were  experienced  with  the  dew  point  meter  which 
occasionally  gave  readings  above  the  air  temperature.  However,  during 
most  measurement  periods  the  dew  point  readings  are  considered  accurate 
and  agreed  very  well  with  measurements  made  at  a  nearby  heliport. 

2.5  Aerosol  Particle  Concentration  and  Size  Distribution 


The  aerosol  concentration  and  size  distribution  were  measured 
with  the  Particle  Measuring  Systems  model  ASASP-300  active  and  model 
CSASP-100-  HV  classical  probes.  This  arrangement  provides  particle 
concentration  measurements  In  eight  overlapping  size  ranges  each 
divided  into  15  resolution  bins.  The  total  range  covers  particles  of 
0.15  pm  to  32  pm  in  diameter.  The  probes  were  installed  on  the  roof  of 
the  receiver  station.  They  were  operated  in  the  standard  axial  fan/ 
accelerator  horn  configuration.  To  prevent  snowflakes  from  being  drawn 
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into  the  sample  cavity,  the  accelerator  horns  were  oriented  downward. 
The  probes  were  calibrated  using  a  Devilbiss  nebulizer  which  produces 
monodispersed  spherical  latex  particles  of  known  diameter. 

2.6  Snow  Characterization 

Samples  of  falling  snow  crystals  were  preserved  for  visual 
inspection  using  the  replica  technique  described  in  Ref.  1.  The  sam¬ 
ples  were  usually  taken  at  half-hour  intervals.  Using  these  samples, 
the  types  of  snow  crystals  were  identified  and  estimates  made  of  the 
mean  particle  radius  T  and  standard  deviation  z.  The  airborne  snow  mass 
concentration  was  also  measured  using  the  Airborne  Snow  Concentration 
Measuring  Equipment  (ASCME)  built  at  DREV  after  the  design  of  Lacombe 
(Ref.  2).  This  device  provides  a  real-time  measurement  of  the  mass 
concentration  of  airborne  snow  which  was  recorded  along  with  the  trans¬ 
mission  measurements.  The  ASCME  is  fully  described  in  Refs.  2  and  3. 

3.0  VISIBILITY  AND  VISIBLE  TRANSMITTANCE 


3 . 1  Observer  Estimates  of  Visibility 


Meteorological  range  or  visibility  is  defined  as  the  greatest 
distance  at  which  it  is  possible  to  see  and  identify  with  the  unaided 
eye:  a)  a  prominent  dark  object  against  the  horizon  sky  during  the  day 

and  b)  an  unfocused,  moderately  intense  light  source  during  the  night 
(Ref.  4).  The  vagueness  of  these  definitions  and  the  subjectiveness  of 
individual  observers  means  that  visibility  estimates  are  generally  not 
very  accurate.  Unfortunately,  in  the  field,  visibility  is  often  the 
only  indication  of  the  degree  of  atmospheric  extinction.  In  the  XSCALE 
module,  visibility  is  an  important  parameter  in  the  prediction  of  IR 
transmittance. 


"N"LAsc. !  - 


An  empirical  formula  that  relates  visibility  V  to  the  vi^l'N 

extinction  coefficient  a  ,  is  Rose hme 1 de r ’ s  formula: 

vl  s 


V 


3.912 

Jvis 


Equation  1  may  be  combined  with  Beer's  law, 
-  La 

t  =  e 


o 


to  obtain  a  relation  between  visibility  and  visible 
i  .e. 


transmittance  t  , 
vis 


V 


3.912  L 


inr 


vis 


[3. 


where  L  is  the  length  of  the  transmission  path.  Figure  3  is  a  scatter 
plot  of  observer  estimates  of  visibility  vs  visibility  derived  from 
visible  transmittance  measurements  using  eq.  3.  The  error  in  the 
observer-estimated  visibility  is  about  352. 


To  eliminate  the  subjective  element  of  observer  visibility 
estimates,  the  XSCALE  snow  algorithm  is  tested  in  Chapter  4.0  using 
values  of  visibility  derived  from  the  measured  visible  transmittance 
using  eq.  3. 

3.2  Extinction  Coefficient  and  Forward  Scattering 

Beer's  law  is  exact  only  for  an  ideal  receiver  with  an  infinite¬ 
ly  small  field  of  view.  For  real  receivers  with  a  finite  field  of  view 
operating  in  a  medium  with  extinction  coefficient  a,  off-axis  contribu¬ 
tions  to  the  detected  light  can  result  in  a  measured  transmittance  that 
is  considerably  greater  than  that  given  by  eq.  2. 
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FIGURE  3  - 


2  3  456789 

CALCULATED  VISIBILITY  (KM) 

Visibility  estimated  by  an  observer  vs  visibility  derived 
from  visible  transmittance  measurements  using  eq.  3 


Snow  particles  are  generally  larger  than  IR  and  visible  wavelengths  and 
thus  a  large  fraction  of  the  energy  scattered  by  snow  is  contained  in  a 
narrow  forward  lobe  (Ref.  5).  Since  the  angular  width  of  the  forward 
scattering  lobe  Is  of  the  same  order  as  typical  transmissometer  receiv¬ 
ers,  forward  scattering  can  contribute  greatly  to  snow  transmittance 
measurements  depending  on  particle  size,  wavelength  and  transmissometer 
design. 
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TABLE  II 

Forward  scattering  correction  to  measured  transmittance 


Measured 

Transmittance 

4  ym 

Corrected  Transmittance 

10  ym  0.633  ym 

1.0 

1.0 

1.0 

1.0 

0.9 

.86 

.88 

.89 

0.7 

.61 

.65 

.68 

0.4 

.29 

.33 

.37 

0.2 

.12 

.14 

.18 

0.1 

.05 

.06 

.07 

The  component  of  measured  transmittance  due  to  forward  scatter¬ 
ing  was  computed  using  the  Mill  and  Shettle  model  (Ref.  6).  The 
"corrected"  transmittance  or  the  transmittance  of  the  unscattered  beam 
was  found  by  subtracting  the  forward  scattering  component  from  the 
measured  transmittance.  Examples  of  measured  transmlttances  and  the 
corresponding  forward  scattering  corrected  values  are  given  in  Table 
II.  The  calculations  were  done  with  an  assumed  snow  particle  radius  of 
0.5  mm.  See  Ref.  1  for  details. 

Table  II  shows  that  the  fraction  of  transmittance  due  to  forward 
scattering  Increases  as  the  transmittance  decreases.  The  average  per¬ 
cent  difference  between  measured  and  corrected  values  is  27Z  at  4  pm 
and  18Z  at  10  ym.  Therefore,  the  contribution  of  forward  scattering  to 
the  measured  transmittance  Is  1.5  times  greater  at  4  ym  than  at  10  ym. 
The  effect  of  forward  scattering  on  the  visible  transmissometer  Is  much 
less  than  that  on  the  IR  transmissometer  with  an  average  percent 
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difference  between  measured  and  corrected  visible  transmittance  of  10%. 
This  has  been  attributed  to  the  fact  that  the  visible  transmissometer 
beam  is  much  smaller  in  diameter  than  the  receiver. 

Extinction  coefficients  calculated  with  eq.  2  using  corrected 
transmittance  measurements  should  be  closer  to  the  true  values  than 
those  derived  from  uncorrected  measurements.  In  Chapter  4.0,  values  of 
visibility  used  by  XSCALE  are  derived  from  corrected  visible  transmit¬ 
tance  measurements. 

4.0  XSCALE  TRANSMITTANCE  PREDICTIONS 


4.1  The  XSCALE  Module 


The  EOSAEL  module  XSCALE  is  a  self-contained  FORTRAN  program 
that  can  estimate  visible  and  IR  transmission  during  low  visibility 
weather  conditions.  XSCALE  calculates  transmittance  through  fog,  rain 
and  snow  using  various  empirical  and  seraierapirical  models.  The  trans¬ 
mittance  component  due  to  background  aerosols  is  calculated  using  Mie 

TABLE  III 


Data  required  for  XSCALE  transmittance  calculation 


Input 

(units) 

Output 

(units) 

Wavelength 

(pm) 

Extinction  coefficient 

(km-1) 

Transmission  path  length 

(km) 

Transmission 

(*) 

Visibi lity 

(km) 

Receiver  radius 

(cm) 

Temperature 

(cC) 

Relative  humidity 

(*) 

* 


H 

V«>1 


M3 

IF 


$**} 

» 
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theory  and  model  aerosol  size  distributions  based  on  size  distributions 
measured  in  many  common  air  mass  types  (Ref.  7).  The  data  required  by 
XSCALE  to  make  horizontal  path  transmittance  measurements  are  listed  in 
Table  III  along  with  the  computation  results. 


The  extinction  coefficient  due  to  the  snow  particles  is 


calculated  using  the  formula 


C 1  *96^  ^  -  0.88k 

a  -  ( - J  (1  +  e  ) 

8  V 


where  k  «  2itr  r  ,/AV  (Ref.  5).  In  eq.  4,  V  is  the  visibility,  r  is 
Q  a 

the  receiver  radius,  X  is  the  wavelength  and  T  is  the  mean  snow 
particle  radius,  r  is  estimated  using  the  empirical  relations  given  in 
Table  IV. 


TABLE  IV 


XSCALE  snow  particle  radius  vs  temperature 


Radius  r  (pm) 

Temperature  T  (°C) 

100 

T  <  -15 

250  +  10T 

-15  <  T  <  0 

250  +  25T 

0  <  T  <  2 

300 

T  >  2° 

% 


>  W  v  • 
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XSCALE  uses  one  of  two  possible  background  aerosol  models  for 
snow  transmittance  estimates.  The  model  used  is  chosen  automatically 
by  XSCALE  based  on  the  input  relative  humidity.  For  humidity  less  than 
95%,  the  aerosol  model  is  one  associated  with  dry  snow  and  contributes 
little  to  visible  or  IR  extinction.  For  humidity  greater  than  95%,  a 
background  of  advective  fog  is  assumed.  The  fraction  of  total  extinc¬ 
tion  due  to  the  fog  depends  on  the  concentration  of  the  fog  particles. 
The  total  extinction  a  may  be  expressed  as  the  sum  of  relative  extinc¬ 
tion  contributions  due  to  fog  and  snow,  i.e. 

a  *  (1  -  B)  a£  +  Ba 
r  s 

where  is  the  fog  extinction  coefficient  and  the  factor  B  is  calcu¬ 
lated  using  the  empirical  relation 

B  »  -  0.025  T  +  0.23  V-1  -  0.021  RH  +  2.47.  [5] 

In  eq.  5 ,  T  is  the  temperature,  V  is  the  visibility  and  RH  is  the 
relative  humidity. 

After  the  extinction  contributions  due  to  snow  particles  and  the 
aerosol  background  are  determined,  the  transmission  along  the  path 
length  L  is  calculated  using  Beer's  law  and  the  calculated  total 
extinction  coefficient. 

4.2  The  Data 


The  visible  and  IR  transmittance  data  used  to  verify  the  XSCALE 
snow  algorithm  were  collected  on  March  1  and  April  2  and  3,  1985.  Time 
plots  of  these  data  are  shown  in  Figs.  4  through  6.  The  extinction  due 
to  molecular  absorption  has  been  removed  from  the  transmittance  data  by 
normalizing  them  to  molecular  transmittances  calculated  using  L0WTRAN 
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as  described  in  Ref.  1.  Thus  the  transmittance  curves  in  Figs.  4-6  are 
due  to  extinction  by  airborne  snow  and  the  aerosol  background.  Figures 
4-6  also  show  the  airborne  snow  mass  concentration  measured  with  the 
ASCME  device  described  in  Refs.  2  and  3.  The  snow  concentration  curves 
indicate  the  intensity  of  the  snow  episode  where  a  snow  concentration 
of  500  mg  per  cubic  meter  of  air  is  considered  a  very  heavy  snowfall. 

The  transmittance  component  due  to  the  background  aerosols  was 
calculated  using  the  measured  aerosol  size  distribution  and  Mie  theory 
as  described  in  Ref.  1.  The  results  show  that  for  the  snow  episodes  of 
April  2  and  3,  attenuation  in  the  visible  due  to  aerosols  is  slight, 
with  the  visible  aerosol  transmittance  being  greater  than  0.97.  IR 
extinction  due  to  aerosols  is  negligible  on  these  days.  The  aerosol 
measurements  indicate  that  on  March  1  there  was  a  higher  aerosol  con¬ 
centration  but  the  transmittance  due  to  background  aerosols  was  greater 
than  0.93  in  all  bands.  The  aerosol  size  distributions  measured 
between  10:00  and  11:00  on  April  2  and  18:00  to  19:00  on  April  3  are 
shown  in  Figs.  7a  and  7b. 

The  temperature,  relative  humidity  and  visibility  data  used  as 
input  to  XSCALE  are  found  in  Table  V  along  with  aerosol  visible  trans¬ 
mittance  x  derived  from  the  aerosol  concentration  measurements, 
aer 

On  March  1  and  April  3,  the  measured  relative  humidity  was  between  99 
and  100%  although  the  values  indicated  in  Table  V  are  95%.  The  reason 
for  this  substitution  was  to  avoid  invoking  the  XSCALE  advective  fog 
background  aerosol  model.  The  aerosol  visible  transmittance  values  in 
Table  V  show  that  the  extinction  due  to  aerosols  was  small  for  all 
three  days.  Therefore,  the  advective  fog  model  which  produces  high 
extinctions  would  be  Inappropriate.  Thus  all  transmittance  predictions 
made  with  XSCALE  use  the  dry  snow  aerosol  model. 

As  discussed  In  Section  3.2,  the  visibility  values  are  derived 
from  the  forward  scatter  corrected  visible  transmittance  measurements 


TRANSMITTANCE  I  «  TRANSMITTANCE 
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FIGURE  6  -  Visible  transmittance,  IR  transmittance  and  snow  mass 
concentration  for  April  3,  1985 

-  Visible - 3-5  pm 

- —8-12  pm  ■  snow  mass  concentration 


4 . 3  XSCALE  Transmittance  Estimates 

Using  the  data  of  Table  V  and  a  path  length  of  0.538  km,  esti¬ 
mates  of  transmittance  at  0.63  pm,  4  pm  and  10  pm  were  made  using 
XSCALE.  For  the  IR  predictions,  the  receiver  radius  was  4.57  cm  and 
for  the  visible  predictions,  the  receiver  radius  was  12.7  cm.  In  Table 
VI,  the  results  are  compared  to  the  transmlttances  measured  at  0.63  pm 
and  in  the  3-5  pm  and  8-12  pm  bands . 

At  0.63  pm  the  XSCALE  transmittance  estimates  agree  very  well 
with  the  measured  values  for  all  three  snow  episodes.  However,  at  4  pm 
the  calculated  transmlttances  are  on  average  151  lower  than  the  meas¬ 
ured  values  for  March  1  and  April  2.  On  April  3,  the  visibility  was 
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TABLE  V 

Meteorological  data  during  snow  episodes 


Time  Visibility 

Temperature 

(°C) 

Relative  Humidity 

(Z) 

T 

aer 

(\“.63  pm) 

March  1,  1985 

10:30 

1.3 

-2.2 

95 

.97 

11:00 

2.7 

-2.1 

95 

.96 

11:30 

.94 

-2.1 

95 

.95 

12:00 

3.1 

-2.0 

95 

.94 

12:30 

2.7 

-2.0 

95 

.94 

13:00 

5.9 

-2.1 

95 

.94 

14:30 

3.2 

-1.8 

95 

.94 

15:00 

1.6 

-1.6 

95 

.94 

15:30 

2.2 

-1.6 

95 

.93 

April  2,  1985 

10:00 

6.4 

-3.3 

86 

_ 

10:30 

6.4 

-3.3 

88 

- 

11:00 

3.1 

-2.7 

86 

.99 

11:30 

5.5 

-2.9 

87 

.99 

12:00 

2.1 

-2.8 

86 

.99 

12:30 

3.3 

-2.8 

85 

.99 

13:00 

6.2 

-2.5 

85 

.99 

Aoril  3,  1985 

17:30 

1.5 

-2.0 

95 

.97 

18:00 

1.0 

-2.4 

95 

.96 

19:45 

1.4 

-2.6 

95 

.97 

20:00 

1.5 

-2.8 

95 

.97 

20:30 

1.1 

-2.6 

95 

- 

20:45 

1.0 

-2.6 

95 

- 

21:00 

1.1 

-2.6 

95 

- 
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TABLE  VI 

Measured  and  calculated  snow  transmittance 


Time 

0.63  pm 

Meas. 

0.63  pm 

Cal. 

3-5  pm 

Meas . 

4  pm 

Cal. 

8-12  pm 

Meas . 

10  pm 

Cal. 

March  1, 

1985 

10:30 

.23 

.24 

.28 

.20 

.23 

.20 

11:00 

.49 

.48 

.48 

.46 

.49 

.46 

11:30 

.13 

.14 

.18 

.11 

.15 

.11 

12:00 

.51 

.50 

.56 

.48 

.57 

.48 

12:30 

.49 

.48 

.56 

.46 

.55 

.46 

13:00 

.72 

.71 

.78 

.70 

.75 

.70 

14 : 30 

.54 

.53 

.59 

.51 

.57 

.51 

15:00 

.30 

.31 

.39 

.27 

.36 

.27 

15:30 

.42 

.41 

.49 

.39 

.49 

.39 

April  2, 

1985 

10:00 

.74 

,73 

.80 

.72 

.77 

.72 

10:30 

.74 

.73 

.84 

.72 

.75 

.72 

11:00 

.51 

.50 

.57 

.48 

.50 

.48 

11:30 

.70 

.69 

.73 

.68 

.73 

.68 

12:00 

.39 

.39 

.48 

.36 

.41 

.36 

12:30 

.56 

.55 

.61 

.53 

.58 

.53 

13:00 

.73 

.72 

.78 

.71 

.77 

.71 

April  3, 

1985 

17:30 

.26 

.27 

.35 

.23 

.30 

.23 

18:00 

.15 

.16 

.28 

.13 

.20 

.13 

19:45 

.24 

.25 

.37 

.21 

.30 

.21 

20:00 

.28 

.29 

.38 

.25 

.33 

.25 

20:30 

.18 

.19 

.30 

.16 

.24 

.16 

20:45 

.14 

.15 

.21 

.12 

.17 

.12 

21:00 

.18 

.19 

.26 

.16 

.21 

.16 
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considerably  less  than  on  the  other  two  days  and  the  estimated  4-pm 
transmittance  is  on  average  41%  less  than  the  measured  values.  The 
results  are  similar  in  the  8-12  pm  band  with  the  average  calculated 
transmittances  less  than  the  measured  values  by  about  20%. 

4.4  Discussion  of  Results 

Atmospheric  transmittance  measurements  made  through  falling  snow 
may  include  a  large  component  due  to  forward  scattered  light.  This 
forward  scattering  component  was  removed  from  the  visible  transmittance 
data  using  the  method  described  in  Section  3.3  of  Ref.  1.  The  visibil¬ 
ity  data  given  in  Table  V  are  derived  from  the  corrected  visible  trans¬ 
mittance  using  eq.  3.  Since  the  XSCALE  estimates  of  visible  transmit¬ 
tance  are  very  close  to  the  measured  values,  it  appears  that  the  XSCALE 
snow  algorithm  has  determined  the  contribution  of  forward  scattering  in 
the  visible  transmittance  measurements. 

TABLE  VII 

Forward  scatter  corrected  (FSC)  transmittance  and  XSCALE  predictions 


Time 

Measured 

4  pm 

FSC 

XSCALE 

Measured 

10  pm 

FSC 

XSCALE 

April  3,  1985 

17:30 

.35 

.24 

m 

CM 

• 

.30 

.24 

.23 

18:00 

.28 

.18 

.13 

.20 

.15 

.13 

19:45 

.37 

.26 

.21 

.30 

.24 

.21 

20:00 

oo 

.27 

.25 

.33 

.26 

.25 

20:30 

.30 

.24 

.16 

.24 

.18 

.16 

20:45 

.21 

.13 

.12 

.17 

CM 

r-M 

• 

.12 

21:00 

.26 

.17 

.16 

.21 

.16 

.16 
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At  4  pm  and  10  pm,  the  XSCALE  estimates  of  transmittance  are 
much  lower  than  the  measured  values.  Table  VII  compares  some  of  the 
XSCALE  estimates  at  4  pm  and  10  pm  with  the  measured  transmittance,  and 
the  measured  transmittance  with  the  forward  scattered  component 
removed . 

At  4  pm,  the  XSCALE  predictions  are  on  average  about  152  less 
than  the  forward  scatter  corrected  values  and  about  402  less  than  the 
measured  values.  At  10  pm,  the  XSCALE  values  are  about  the  same  as  the 
forward  scatter  corrected  values  and  23%  lower  than  the  measurements. 
Since  the  XSCALE  transmittance  predictions  are  consistently  less  than 
the  measured  values  but  close  to  the  forward  scatter  corrected  values, 
it  appears  that  the  XSCALE  snow  algorithm  underestimates  the  contribu¬ 
tion  of  forward  scattering  in  the  IR  transmittance  measurements. 

According  to  Table  IV,  which  gives  7  as  a  function  of  tempera¬ 
ture,  the  mean  snow  particle  radius  used  by  XSCALE  in  the  calculations 
was  about  230  pm  for  all  three  snow  episodes.  Values  of  snow  particle 
radius  determined  by  analyzing  replicas  of  the  snow  crystals  taken 
during  each  snowfall  are  given  in  Table  VIII  with  the  standard 
deviation. 

TABLE  VIII 

Observed  snow  particle  radius  7 


Episode 


r  ( pm) 


Std.  dev. 


March 

1. 

1985 

1600 

1100 

April 

2, 

1985 

900 

700 

April 

3, 

1985 

800 

500 
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A  meaningful  representative  radius  for  snow  particles  is  diffi¬ 
cult  to  ascertain  since  the  range  of  particle  sizes  is  very  large  and 
the  particles  are  irregularly  shaped.  However,  Table  VIII  shows  that 
the  mean  radius  of  snow  particles  observed  during  the  three  episodes 
are  much  larger  than  the  230  pm  used  in  the  XSCALE  calculation.  Fur¬ 
thermore,  when  using  the  Mill  and  Shettle  forward  scattering  calcula¬ 
tion,  as  described  in  Ref.  1,  optimal  results  were  achieved  using  r  =* 
500  pm.  Use  of  an  uncharacteristically  small  value  for  r  by  XSCALE 
could  explain  why  the  algorithm  underestimates  the  forward  scattering 
contribution  to  the  measured  transmittance. 

5.0  THE  XSCALE  SNOW  ALGORITHM 


5.1  Effective  Extinction  Coefficient 


An  ideal  transmissometer  has  a  receiver  with  an  infinitely  small 
field  of  view.  Transmittance  measurements  made  with  such  a  transmis¬ 
someter  are  related  to  the  volume  extinction  coefficient  of  the  medium 
exactly  according  to  Beer's  law, 


where  a  is  the  extinction  coefficient  and  L  the  transmission  path 
length.  For  a  receiver  with  a  finite  field  of  view,  the  measured 
transmittance  may  be  larger  due  to  off-axis  light  scattered  into  the 
receiver.  The  extinction  coefficient  calculated  using  such  measure¬ 
ments  and  eq.  6  is  thus  the  apparent  or  effective  extinction  coeffi¬ 
cient  a,  which  is  smaller  than  a  depending  on  the  fraction  of  forward 
scattered  energy  collected  by  the  receiver. 

In  order  to  estimate  transmittance  through  falling  snow,  the 
XSCALE  model  calculates  the  effective  extinction  coefficient  a  and 
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derives  the  transmittance  from  eq.  6.  A  derivation  of  the  XSCALE  snow 
algorithm  developed  by  Seagraves  (Ref.  5)  is  summarized  in  the  follow¬ 
ing  section. 

5.2  Derivation  of  the  Snow  Algorithm 

The  detection  of  forward  scattered  light  by  a  receiver  contrib¬ 
utes  to  the  measured  transmittance  and  makes  the  particles  appear  to  be 
less  effective  scatterers  or  to  have  a  smaller  extinction  efficiency. 
Seagraves  defines  an  effective  extinction  efficiency 

Qef f  “  Q  -  AQ  [7] 

where  AQ  is  the  reduction  in  Qe^£  due  to  forward  scattering.  As  shown 
in  Fig.  8,  a  receiver  of  radius  r^  subtends  an  angle  0q  at  a  point  i  on 
the  beam  axis. 

For  a  particle  at  i,  AQ  is  found  by  integrating  the  phase  func¬ 
tion  p(0)  of  the  particle  over  the  solid  angle  subtended  by  the 
receiver,  i.e. 


0 

o 

AQ(0q)  *2 it  /  p(0)  sln9d0 .  [8] 

o 

The  phase  function  used  to  evaluate  eq.  8  is 

2  2J1(nsino)  2 

p(9,n)  - —  ( - j 

4tt  nsin0 

where  n  is  the  particle  size  parameter  equal  to  2rr/X  and  J1  is  the 
first-order  Bessel  function  of  the  first  kind.  p(0,n)  is  an  approxi¬ 
mate  phase  function  derived  from  Hie  theory  and  applies  to  homogeneous 
spheres  larger  than  the  wavelength  (Ref.  8). 


9 
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•FIGURE  8  -  Geometry  used  to  calculate  forward  scattering  contribution 
in  XSCALE  snow  algorithm 


Hodkinson  found  that  the  forward  diffraction  lobe  for  randomly 
oriented  nonspherical  particles  was  similar  to  that  of  spheres  of  equal 
geometric  cross  section  (Ref.  9),  so  that  p(9,n)  may  be  considered  a 
reasonable  approximation  to  the  phase  function  of  snow  particles. 
Evaluating  the  integral  in  eq.  8  yields 

2 

AQ(9q)  ■  1  -  j2(nsin9Q)  -  JjCnsinB^.  [  1 0  j 

Since  sin9  “  r,/£AQ(9  )  is  a  function  of  distance  along  the  transmis- 
o  a  o 

sion  path  and  receiver  radius,  forward  scattering  contributions  from 
particles  distributed  along  the  beam  axis  may  be  accounted  for  by  cal¬ 
culating  a  mean  extinction  efficiency  This  is  found  by  integrat¬ 

ing  AQ  with  respect  to  £  along  the  beam  path  and  subtracting  it  from 
Q(A,r),  i.e. 

,  L 

\{{  m  QU.r)  -  r  /  AQ(£)d£,  [11  ; 

o 
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where  L  is  the  length  of  the  transmission  path.  Since  snow  particles 
are  generally  larger  than  visible  and  IR  wavelengths,  Q(A,r)  is  taken 
to  be  equal  to  2.  Seagraves  gives  an  approximate  solution  to  eq.  11 
as 


nr 


-  0.88 


‘eff 


=  e 


) 


+  1. 


[12] 


Note  that  as  r^  becomes  vanishingly  small,  TJ  ^  becomes  equal  to 
2  which  means  that  this  model  predicts  that  forward  scattering  effects 
are  negligible  for  a  sufficiently  small  receiver.  The  effective 

k  _ 

extinction  coefficient  a  is  then  found  by  substituting  Q  ,,  in  the 

eff 

definition  of  the  extinction  coefficient. 


a 


Q  I 


dN 

IF 


itr2dr. 


where  N  is  the  number  density  function  of  the  particles.  Combining 
eqs .  12  and  13  yields 

r  ("  0.88  k)  -I  r  dN  2j  r  i  /  1 

o  -  [e  +  !]  J  *r2dr  [14] 

o 

where  k  ■  nr^/L.  Using  eq.  14,  the  ratio  of  a  for  any  two  wavelengths 
may  be  written 


<*1 


1  +  e 


(-  0.88  kx) 


<*2  j  +  k2^ 


15 


In  order  to  relate  a(A)  to  visibility,  the  of  the  human  eye 
is  taken  to  be  equal  to  2  since  it  Is  a  very  small  receiver.  The  visi¬ 
ble  extinction  coefficient  is  given  by  Koschmeider' s  formula  as 


■I 


1 


I 


» 

3 

t 

t- 

1 


b 


J 

i 

■ 

m 

3 

k> 
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3.912 


vis 


[16] 


Combining  eqs.  15  and  16  yields  the  effective  extinction  coefficient 


1.96  .  -  0.88k. 

a(k,V)  =  -  (1  +  e  ). 


[17; 


Substituting  a  in  Beer's  law  we  obtain  the  XSCALE  snow  transmittance 
algorithm 


T(k,L,V)  =  exp  {-  —  (1  +  e 

V 


0.88k 


)}• 


[18] 


In  the  XSCALE  code  (EOSAEL  version  1984),  the  coefficient  k  has  the 

value  nr  ,/V. 
d 

5.3  Detector  Diameter  vs  Beam  Diameter 


The  effective  extinction  efficiency  defined  in  eq.  11 

accounts  for  reduced  extinction  due  to  light  forward  scattered  within  a 

solid  angle  tt  r2  /l2  about  the  forward  direction  as  shown  in  Fig.  8. 
d 

This  means  that  only  contributions  from  scatterers  near  the  beam  axis 
are  accounted  for.  For  a  transmissometer  whose  beam  width  is  smaller 
than  the  diameter  of  the  receiver  this  should  be  a  realistic  model.  In 
Chapter  4.0,  XSCALE  visible  transmittance  estimates  are  shown  to  agree 
very  well  with  measured  values.  The  visible  transmittance  was  measured 
using  a  narrow  beam  transmissometer  with  the  beam  diameter  one  tenth 
the  receiver  diameter.  The  good  agreement  between  predicted  and  meas¬ 
ured  transmittance  in  the  visible  indicates  that  the  model  is  a  good 
representation  of  a  narrow  beam  transmissometer. 

As  discussed  in  Section  3.2,  the  visibility  data  used  for  the 
XSCALE  calculations  were  derived  from  the  measured  visible  transmit¬ 
tance  from  which  the  forward  scattered  contribution  was  removed  using 
the  method  of  Mill  and  Shettle  (Ref.  6).  The  fact  that  the  XSCALE 
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visible  transmittance  predictions  are  very  close  to  the  measured 
visible  transmittance  means  that  the  XSCALE  estimate  of  the  forward 
scattering  contribution  in  the  visible  agrees  very  closely  with  the 
Mill  and  Shettle  model. 

For  transmissometers  with  a  beam  much  larger  than  the  receiver, 
forward-scattered  energy  may  be  received  from  particles  far  from  the 
beam  axis  as  illustrated  in  Fig.  9.  It  seems  likely  that  for  wide-beara 
transmissometers  the  effective  extinction  efficiency  calculated  in  eq. 
11  will  be  too  large  since  only  the  forward  scattering  contributions  of 
particles  near  the  beam  axis  are  acknowledged.  This  is  supported  by 
the  fact  that  the  XSCALE  predictions  in  the  IR  were  consistently  below 
the  measured  values  in  both  the  3-5  pm  and  8-12  pm  bands,  suggesting 
that  XSCALE  underestimated  the  magnitude  of  the  forward-scattered 
transmittance  component. 


FIGURE  9  -  Transmissometer  with  beam  much  larger  than  the  receiver 
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6.0  CONCLUSIONS 


Estimates  of  transmittance  through  falling  snow  made  by  the 
EOSAEL  module  XSCALE  were  compared  to  measurements  of  visible  and  1R 
transmittance  made  at  DREV.  The  XSCALE  algorithm  requires  as  input  the 
visibility  which  was  derived  from  measurements  of  visible  transmittance 
using  the  Beer  and  Koschmieder  laws.  The  visible  transmittance  values 
used  to  calculate  the  visibility  were  corrected  for  forward  scattering 
effects  using  the  Mill  and  Shettle  algorithm  (Ref.  6).  Using  the 
derived  visibility  values  as  input,  the  XSCALE  predictions  of  visible 
snow  transmittance  agreed  very  well  with  measurements  made  using  a 
narrow-beam  laser  transmissometer  with  a  receiver  larger  than  the  beam. 

XSCALE  predictions  of  snow  transmittance  in  the  IR  were  on 
average  about  30%  less  than  the  measured  values  in  the  3-5  pm  band  and 
about  20%  less  than  measured  values  in  the  8-12  pm  band.  This  could  be 
due  in  part  to  the  fact  that  the  XSCALE  algorithm  uses  a  mean  particle 
radius  which  is  much  smaller  than  those  observed  during  the  snow 
transmittance  measurements.  The  IR  transmittance  measurements  were 
made  using  a  transmissometer  with  a  divergent  beam  much  wider  than  the 
receiver  at  the  receiver  end. 

It  has  also  been  shown  that  in  calculating  the  forward-scattered 
transmittance  component,  the  XSCALE  algorithm  does  not  take  into 
account  the  contribution  of  scatterers  found  far  from  the  beam  axis. 
This  could  explain  why  XSCALE  underestimates  the  forward-scattered 
component  for  transmissometers  with  a  beam  much  wider  than  the  receiver 
since  scatterers  far  from  the  beam  axis  can  contribute  greatly  to  the 
energy  collected  by  the  receiver.  On  the  other  hand,  the  fact  that  the 
XSCALE  predictions  are  very  close  to  the  measured  visible  transmit- 
tances  shows  that  the  algorithm  is  a  good  representation  of  a  narrow- 
beam  transmissometer  whose  beam  is  much  smaller  than  the  receiver 
diameter. 
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